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The exact fermion master equation previously obtained in [Phys. Rev. B 78 , 235311 (2008); 
New J. Phys. 12 , 083013 (2010)] describes the dynamics of quantum states of a principal system 
of fermionic particles under the influences of external fermion reservoirs (e.g. nanoelectronic sys¬ 
tems). Here, we present the general solution to this exact fermion master equation. The solution 
is analytically expressed in the most intuitive particle number representation. It is applicable to 
an arbitrary number of orbitals in the principal system prepared at arbitrary initial states. We 
demonstrate the usefulness of such general solution with the transient dynamics of nanostructured 
artificial molecules. We show how various initial states can lead to distinct transient dynamics, 
manifesting a multitude of underlying transition pathways. 
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I. INTRODUCTION 

Nanoelectronic systems are versatile platforms to ex¬ 
plore quantum device applications. Conventionally, the 
Schwinger-Keldysh nonequilibrium formalism is applied 
to study their quantum transport propertiesji where the 
quantum states of the devices are not necessarily tackled. 
On the other hand, for the interests of quantum informa¬ 
tion processing, directly addressing the quantum states 
of device systems has been fundamentally important)^ 
Many experimental endeavors have gone to tailor and 
manipulate quantum states of the targeted devicesi^ 
Quantum devices are open systems where interactions 
with environments can significantly affect the states of 
the devicei^^— Equations of motion that govern the time 
evolutions of the states of open systems, known as mas¬ 
ter equations, then become common theoretical tools. To 
have complete characterization of the device operations, 
solutions to master equations, namely, explicit expres¬ 
sions of the quantum states in terms of the reduced den¬ 
sity operators, are strongly sought for4“— 

Preparing as well as reading quantum states have 
been experimentally demonstrated with atoms, ions and 
moleculeSfi^^— superconducting circuits^^— photonic 
modes^^— and electrons in quantum dots (QDs)^^— 
just to name a few. Reconstructing quantum states 
from experimental data of measurements of observables, 
known as quantum state tomography, have also been 
intensively investigated on the physical and algorith¬ 
mic aspects4^“— On the intersections of quantum in¬ 
formation processing, nonequilibrium phenomena and 
nanoscale electronics lies a class of nanoelectronic sys¬ 
tems. They consist of a central system containing a num¬ 
ber of discrete states, exchanging electrons with electron 
reservoirs, as electrodes or other continua in host ma¬ 
terials. Sophisticated nano-fabrication technology makes 
these systems plainly available and highly tunable. Alter¬ 
natively, using cold atoms with optical controls, fermionic 
atoms with discretized montional stateas in a trap, con¬ 
necting to macroscopic fermion reservoirs, have been real¬ 
ized experimentally as analogs to the aforementioned na¬ 


noelectronic systemsi^^— Their atomic transport prop¬ 
erties have also raised much theoretical interests 
Therefore, exploiting the quantum states of these sys¬ 
tems for both purposes of device-oriented and fundamen¬ 
tal studies is very attractive. 

In addition to the Schwinger-Keldysh formalism, mas¬ 
ter equation approaches have also been widely ap¬ 
plied in nanoelectronic device systems as well.— 
Perturbations on tunnelings to leads are often em¬ 
ployed when Coulomb interactions within the cen¬ 
tral region are considered— For effectively non¬ 
interacting fermions, exact master equations have been 
obtaine d^^d^ via the Feynman-Vernon’s influence func¬ 
tional approach^ and strong Coulomb repulsions can be 
taken into account by excluding double occupancies via 
modifying the rate equationsJ^ It is known that the co¬ 
efficients in the master equations reveals information of 
decoherence properties— However, to have direct 
access to the quantum states, solutions to the master 
equations are demanded. Solving the master equations 
in real time for the full non-Markov dynamics is gener¬ 
ally challenging. Here we present the general analytical 
solution to the exact master equations for effectively non¬ 
interacting fermions. We provide direct and complete de¬ 
scriptions of the quantum states in terms of the reduced 
density operator in the Fock basis. 

As a demonstration, we apply this general solution to 
the study of transient dynamics of nanoscale artificial 
molecules. This could be a double-well trap for fermionic 
atoms or a coupled double-quantum-dot (DQD). For such 
two-orbital systems, we have previously obtained the an¬ 
alytical solution only under the condition that the initial 
state is an empty statei^ Here we have arrived at analyti¬ 
cal solutions that enable us to explore different transition 
pathways arising from different initial states. 

The organization of this article is the following. In 
Sec. Ill Al we first briefly review the exact master equation 
for a class of nanoelectronic systems and their alternative 
realizations using cold fermion atoms. Our main result, 
the analytical expression of the reduced density matrix 
that specifies the quantum state of the principal system, 
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is prescribed in Sec. imi In Sec. nni we illustrate the 
utilization of such solution with the transient dynamics 
of artificial molecules. A summary is made in Sec. IIVI 

II. EXACT FERMION MASTER EQUATIONS 
AND THEIR GENERAL ANALYTICAL 
SOLUTION 

A. Exact fermion master equations 

We consider the following general Hamiltonian for 
fermions, 

^.{t) ='Hs{t) (la) 

where 

^ £ij{t)a\aj, (lb) 

z.jGS 


reservoirs, each with its own initial chemical potential 
and temperature T^, namely, 

Ptot (io) = /5 (to) ® n 

CK 

where 

- (f exp [- jn^ito) - /ksTg] 

tr exp [- {HM - I^bTo\ ’ ^ ^ 

with Ha(t) = Efceaeafe(t)cifeC„fc, = _ Efcea , 

and fcs, the Boltzmann constant. The initial state of the 
principal system, /3(to), is arbitrary. 

The exact equation of motion describing the time evo¬ 
lution of the reduced density operator is readily given 
by^^ 


is the Hamiltonian for the principal system, and 

ttE(t) = X! 51 (Ic) 

Oi k£a 


is the sum of Hamiltonians of fermion reservoirs, each la¬ 
beled by a. The exchange of particles between the prin¬ 
cipal system and the reservoirs is described by 

'^T(t) = 55 + yic^k{t)alc^^'^ (Id) 

i^ock 


Here the subscript iGS = {l,2,---, D} enumerates the 
orbitals (discrete states) within the principal system and 
k G a stands for the continuum states within that reser¬ 
voir a. The fermion field operator ai(a|) or Cq,j,(c)j,j,) then 
annihilates (creates) a particle on orbital i € S or state 
k G a. The hopping amplitude between an orbital i in the 
principal system and a state k in reservoir a is given by 
Vaki{t) = Viaki^)- The indices i and k include both exter¬ 
nal (motional) as well as intrinsic (spin or pseudo-spin) 
degrees of freedom of the fermion in question. The time- 
dependence of the Hamiltonian parameters can arise from 
externally applied fields (via gate-voltages in the context 
of nanoelectronics or laser fields in the context of their 
cold-atom analogs). 

The quantum state of the principal system as an open 
system is completely given by the reduced density oper¬ 
ator, defined as 


p (t) = trEPtot (t) , (2) 

where ptot (t) is the density operator of the whole system 
at time t and trE means tracing over all reservoir degrees 
of freedom. 

We assume as usuat^ i^^’^^i^^ that the total system at 
initial time t = to is in a state ptotito) which is a dis¬ 
entangled product of thermal equilibria of the separate 


where the superoperators C^{t) are expressed explicitly 

by 


= - 55 [aUiKt) + 4p(t)o-3] 

+ B,aij {t)ajaj pit) +h.c.'^, 

C(i)/5(0 = 55 [ajKt)4 + P(<)«j«l] 

'ij 

+ Kmj (t)ajP(t)a|-I-h.c.j, (5) 

The time-dependent dissipation and fluctuation coeffi¬ 
cients in Eqs. (IS|), Ka(t) and Xait), are explicitly deter¬ 
mined by the elementary functions, it(r, to), tt(u0 
t;(r, t), via the relations 

Ka{t)= dTgait,T)u{T,to)u~^it,to), (6a) 

Jto 

Xait) = / dr {gait,T)viT,t) -gait,T)uiT,t)} 

Jto 

— K,ait)vit,t). (6b) 


These elementary functions obey the following 
dissipation-fluctuation integrodifferential equations 
of motion 

d r 

—u (r, s) + iE (t) u (t, s)+ dr'g (r, r') u (t', s) = 0, 

(7a) 

for to < s < T < t, 

t;(T, t)= f ds f ds'tt (r, s) g (s, s') it (s', t), (7b) 

J to J to 

and m(t, t) = [^(t, r)]f with m(s, s) = Id- They are 
directly related to the nonequililbrium Green functionsJ^ 
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Here [E {T)]ij =e.^.(r) and 

9{tl,t2) ='^ga{tl,t 2 ), g{tl,t 2 ) = '^ga{tl,t 2 ) (8a) 

a a. 

with 

ga{ti,t 2 ) = J (8b) 

ga{ti,t2) = y ^/a(w)ra (w,tl,t2)e““^*l“*=\(8c) 

are the reservoir correlation functions and fa{i^) is the 
fermi distribution function for the initial equilibrium of 
the reservoir a. The time-dependent spectral density 
(level-broadening) function is given by 

[Tq {U},ti,t2)]ij = 

Viak ih) (t2)S (w - Sakito)) ■ 

kGcx 

(8d) 

B. General analytical solution 

The above master equation was derived using the 
Feynman-Vernon’s influence functional approach in the 
fermion-coherent-state representation. A matrix element 
of the reduced density operator in the fermion-coherent- 
state representation is expressed a a^^d^ 

(^/|p(OII/)^d/x(|o)dAi(^o )/(!/, I/, t||o,ioX^o|p(io)|lo), 

(9a) 

where the fermion coherent states are given by aj^) = 

CilC) and (||a| = with | = (fi,C 2 ,--- ,^Df be¬ 
ing a column vector of Grassman numbers and 

I = being a row vector. The propagat¬ 

ing function in the fermion-coherent-state representation 
is given by 

J V (9b) 

in which 5c [^,4] is the action of the principal system 
and is the influence functional obtained by 

integrating out the reservoirs’ degrees of freedom (see 
Refs. [I^fl3| for the explicit derivations). Specifying the 
action of the principal system with Eq. m. the Grass- 
mann number integral in Eq. Ilbl has been exactly car¬ 
ried out, yielding, 

J (if, I/, i|lo, io) = det [Id - v (t)] x 

exp[^fJfo{t)^o+ioJof{t)^'f+ifJffit)^'f+ioJoo{t)^o) , 

(9c) 


in which we have introduced D x D matrices, 

Jfo (t) = [Jo/ (t)]^ = (Id - v{t))-^u{t), 

Joo (t) = M'l'(t)(lD - v{t))~^u{t) - Id, 

Jff {t) = i'i-D - - Id, (9d) 

with u{t) = u{t,to) and v{t) = Since fermion co¬ 

herent states are used only as a mathematical tool, it is 
not straightforward to gain understanding of the physical 
properties of the resulting quantum states p{t) directly 
from Eq. (I9|). To meaningfully express the reduced den¬ 
sity operator, we have to work in the Fock-state repre¬ 
sentation. 

The Fock basis is defined by 

,nD) = (a|)"i(4)”" •••(a|,)”^|0), (10) 

where |0) is the conhguration where the principal sys¬ 
tem is completely unoccupied. Here ^ {Oj 1} 

enumerates the number of particles occupying the fth or¬ 
bital in the principal system, for i € S. However, when 
D is arbitrary, it is not convenient to denote these states 
using a set of D integers as done in Eq. cni). An n- 
particle configuration is specified by a certain set of n oc¬ 
cupied orbitals, = {ii, - ■ ■ , i„}, in which ii, • • • , € 

S. Ambiguity of defining its corresponding n-particle 
state arises from the anti-commutating property of the 
fermion field operators, namely, • • • a|^a|^a|j0) = 
—|0)- To avoid this ambiguity, we enforce 
an ordering regulation to all many-particle states using 
the notation, 

where P stands for an orderer that = 

, Pi^\ is an ordered set, in which all 

the entries Pil^\--- ,Pi^\Pi^j^^ £ are with the 
ordering Pi^^ > ••• > > Pi^i \ With Eq. (fTTI) . 

we then unambiguously write down the elements of the 
reduced density matrix as 

(Pi(-)|p(t)|Pj(")), (12) 

to represent the coherence between the configuration 
and if ^ or the probability for the system 
to be in the conhguration if = j^'^\ 

Using the transformation between the Eock states and 
the fermion coherent states^ 

(Pi(”)||) = (13) 

Eq. ([H]) can be cast into the Eock-state representation for 
arbitrary D and arbitrary initial state of the principal 
system. This reads, 

(Pi(")|p(t)|Pr(”)) = det [Id - v (t)] 

VUE [jU’;i"f,(<)(f’«'’"»ip(i»)|p!>'’"V. 

m=0 Pa(^) Pbi^) 

(14a) 


X 
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where det [Ld — v (f)] > 7 ^(m )(0 connects the initial 
value of the reduced density matrix evaluated between 
the configurations and 6 ^™^ to the later value of the 
reduced density matrix evaluated between the configura¬ 
tions and r^'^\ Explicitly, they read 


j{n^Q) (n^Q) / r -r / 

'7g(o) g(d) (Z) = det (JJff (t)]j(n) J , 

= (-1)’” det ([Joo (t)]£,’) , 




(t) = (-l)™det 


[•^00 (i)]b(m) 

[jfb (oirdT' 


[Jo/ miZ] 

(14b) 


where represents the vacuum state. The notations 

[Joo (Z)]£,\ [Jo/ (Z)];':,, [J/O (OirdT’ and [Jff (/)]£’ in 
Eq. (Il4bl) are m x m, m x n, n x m and n x n matrices 
respectively, whose elements are given by 


([Joo (l)]b(m) 

) 

0 

II 

n)p„(™ 

([Jo/ miZ 

) 

— [Jo/ it)] 

")p^(": 


/ tj 


3 

([j/o it)]:L 7 

), 

a 

0 

II 

•) p„(m: 


/ tj 


3 

{[jffit)C' 

) ij 

= [J// 

.)p^(n) 

3 


(14c) 


for i,j being in their respective ranges. 

Equation (11^1) is the general analytical solution to the 
exact fermion master equation. Note that due to the 
conservation of the total particle number, the dynam¬ 
ics for certain elements of the reduced density operator, 
|;o(t)|P 6 *^™i) with n m is decoupled from those 
with n = m. Since there is no pair creation or pair 
annihilation present in the Hamiltonian of Eq. , we 
are not interested in these matrix elements resulted from 


J 


superposition between configurations of different parti¬ 
cle numbers. For physically meaningful preparations of 
the initial states, the elements {Pl^'^'^\p{to)\Pb^'^'>) with 
n ^ m are zero for all configurations and 6 ^™^ and 
all values of n and m. The subsequent time evolution 
therefore guarantees (PZ("’)|;o(t > to)\Pb^"^'>) = 0 for any 
Z*^") and Z)^™) with n ^ m. 


III. EXAMPLE: TRANSIENT DYNAMICS OF 
ARTIFICIAL MOLECULES 


In this section, we demonstrate the convenience of us¬ 
ing the above solution with an example, namely, the 
dynamics of an artificial molecule formed in electronic 
nanostructures or its atomic analog. It consists of two 
atomic orbitals (localized orbital in a QD or a potential 
well) tunnel coupled to each other. This corresponds to 
set = 2 in Eq. CEI. We denote the possible configura¬ 
tion of occupying the atomic orbitals by | 0 ) (both atomic 
orbitals are not occupied), | 1 ) (orbital 1 is occupied and 
orbital 2 is empty), | 2 ) (orbital 2 is occupied and orbital 
1 is empty), and |3) (both atomic orbitals are occupied). 

The expression for the reduced density operator in the 
Fock basis, Eq. m, then explicitly reads, 

(*l/3(i)|j)=p..w=pifw+prw> (15a) 

for i,j € {0,1,2,3}. The part that does not depend on 
the initial state of the principal system is 

plTit)=det [v (Z)], 

PiTit) = vii (Z) - det [v (Z)], 

P27'{t) = V 22 {t) - det [i; (Z)], 

Po”'(0 = 1 - (i'll (.t) + V 22 {t)) + det [i; (Z)], 

PlT'it) =V2i{t). (15b) 

The part that contains the initial state is given by 


Pa”*' (t) =Vll (t) g °22 (t) + V 22 (t) Pii (z) - Vi 2 (Z) (Z) - U 21 (Z) (0 

+ P 33 {h) det [m'I' (Z) u (Z)] , 

PlT*' it) = (1 - 'V22 (t)) Qii (Z) - Vii (Z) 0^2 (t) + Vl2 (t) 021 (t) + V21 (t) 0?2 (t) 

- P 33 (to) det [m'I' (Z) u (Z)] , 

P22‘*' it) = (1 - it)) 022 it) - V22 it) 0?i (Z) -b i;i2 (Z) 0^i (Z) -f 1121 (Z) 0?2 (t) 

- P 33 ito) det [«■!■ (Z) u (Z)] , 

Pto*' it) = - [(1 - V22 it)) 0?1 (Z) + (1 - 1^11 it)) 0^2 it) + Vl2 it) 021 it) + V21 (Z) 0?2 (*)] 

+ P 33 ito) det [ut (Z) u (Z)] , 

Pm*' it) =021 it) , (15c) 


where the dependencies on the initial state of the molecule appear explicitly with 033 (Zq) and 

4 it) = X! [“ it)]iMaip{to)al,] [u (Z)]!,^. . (15d) 

Li'eS 
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We have denoted Vij{t) = [v {t)]ij for i,j G {1,2}. In 
case that the initial state of the principal system is an 
empty state, then Eq. CSl) beetles (t) = (t), 

reproducing the solution in Ref. 

In what follows, we concentrate on the role of the initial 
states on the dynamics of the molecule. We consider 
the situation where the total Hamiltonian is independent 
of time. Then the spectral density in Eq. ([S]) becomes 
independent of time, namely, 

[Ta (w,tl,t 2 )]ij = r“(w) = 2t:^5(uJ - eak)ViakVakj- 

k^OL 

(16) 

We use the widely applied wide-band limit, leading 
Eq. ([HI) to r“ (w) = Eg. In this case, Ihut^oo u (t) = 0, 
and we immediately see from Eq. (USD that the part 
which carries the information of the initial states of the 
molecule, vanishes in the steady limit 

hm pjf ■ (t) = 0. (17) 

The initial state of the molecule only has a transient ef¬ 
fect. We also assume there is only one reservoir so that 
the label a in Eq. (USD can be removed. 

The molcular bonding and anti-bonding states (BS and 
AS) are obtained with equal on-site energies eu = £22 = 
Eq. In this basis, the Hamiltonian of the molecule be¬ 
comes 

Tis = H(d^gdAs ~ c^jas^Bs)! (18a) 

where the molecular orbitals in terms of atomic orbitals 
are 

^BS “ 

4 s = -^(«i-4)> (18b) 

and fl = |£ 2 i| is the bonding strength between the two 
artificial atoms. The probabilities {B\p{t)\B) = Pgg{t) 
and {A\p{t)\A) = p^^(t), where |i?) = dgslO) and |A) = 
d^g|0), tell to which extent the principal system is in a 
well-defined molecular state. 

There are two categories of transition path ways that 
lead the changes of the molecular quantum states. One is 
the internal transition induced by coherent coupling be¬ 
tween the two atomic orbitals. The other is via exchang¬ 
ing particles with the external reservoir. By preparing 
the molecule at different initial states, we are able to un¬ 
fold the actions of these different transition pathways. 

In Fig. |T] we respectively monitor the time evolution 
of probabilities for the one-particle BS and AS, starting 
from various initial states of the molecule. The situation 
where BS and AS are well separated in energy, namely, 
n > 2r, where E = (Eu -|- E 22)/2 is the average level 
broadening, is shown on the left panel. The situation 
that their separation is within the broadening B < E is 


displayed on the right panel. Since BS and AS are the 
eigenstates of the bare molecule, the transition |A) —|i?) 
occur only via exchanging particles with the reservoir, 
through |A) —>• |0) —>■ \B) and |A) —^ |3) —|i?). This 
is true also for the transition \B) —>■ |A). Therefore the 
resulting probability Pgg{t) {Paa^^)) obtained by start¬ 
ing the system from the state |A) (|i?)) in the transient 
process is not as high as those obtained by starting the 
system from |0) or |3). This is shown in Fig. [U where 
the value of the red solid curves (for pBB{t) starting 
from AS in Fig. [TDaI),(a2) and for pAA{t) starting from 
BS in Fig. [TDbI),(b2)) are transiently smaller than the 
other curves (for those starting from empty state or two- 
particle state). When AS and BS are well separated in 
energy with the chemical potential of the reservoir be¬ 
ing placed below AS but above BS, then the molecule 
tends to reside in the BS, leaving negligible value of 
PAA(t). Such preference is shown in Fig. UDal) (with 
larger pBB{t)) and in Fig. |TDa2) (with smaller pBB{t))- 
The probability for AS can be enhanced by reducing B to 
be smaller than the broadenings (see Paa(^) Fig.|TDb2) 
being larger than those in Fig. jTDbl)). 


The probabilities for empty state and two-particle 
state, starting from other initial states, are shown in 
Fig. 12 The transitions |0) —>■ |3) and |3) —>• |0) involve 
intermediate one-particle states, |0) —>■ \B{A)) -A |3) and 
|3) —>■ \B{A)) -A |0). Therefore it is less efficient to reach 
jo) (1 3)) from |3) (|0)) than to reach it by starting the 
system from the one-particle states \B) and |A). This 
is shown in Fig. |2 where the values of red solid curves 
are transiently smaller than that of the others. Since the 
one-particle BS is less preferable in case of B < E than 
the case with B > 2E, the probability to stay out of one- 
particle space is relatively higher in the former case (see 
Pooit) ( 7 * 33 ( 0 ) iFe right panel is greater than that in 
the left panel in Fig. |2for later times). 


The different effects due to external reservoir and due 
to internal coherent coupling between the atoms can be 
made more distinct. We monitor the dynamics of prob¬ 
abilities 7 * 22(0 s-iib Pii(0 foi' the localized one-particle 
states in Fig. 12 In Fig. EDa) we show the time evolu¬ 
tion of 7 * 22(0 resulting from initializing the molecule in 
states |0) and |3). In contrast to these smooth evolution 
processes in Fig. EDa), clear oscillation of P 22 (t) in time 
is observed if one starts with one particle localized in |I) 
(see the red solid curve in Fig. Elb)). The oscillation of 
Pii{t) starting from the same state is shown to be com¬ 
plementary to the oscillation of P 22 (0 (see the blue solid 
line marked with white circles in Fig. EDb)). This mani¬ 
fests the coherent tunneling between the two atoms. We 
thus demonstrated the convenience of directly applying 
Eq. (1141) to study the transient dynamics resulted from 
various initial states. 
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FIG. 1: (color online) The left panel (plots (al) and (bl)) 
is for fl = Sr and the right panel (plots (a2) and (b2)) is 
for fl = r/5. Plots (al) and (a2) show how Pbb{^) evolves 
in time starting various initial states as labeled. Plots (bl) 
and (b2) show the time evolutions of Paa(^)- The parameters 
used here and also in other figures are p. = eo and fcsT = O.IP 
with r 22 /rii = 0.7 and ri 2 = 0.2y/TuT^. 
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FIG. 2: (color online) Plots (al) and (a2) show the time evo¬ 
lutions of Poo (t) initiated from various states as labeled. Plots 
(bl) and (b2) display the time evolutions of Pssit). In plots 
(al) and (bl) we let Q = SP. In In plots (a2) and (b2) we let 
fi = r /5 


IV. SUMMARY 


FIG. 3: (color online) Plot (a) presents the time evolutions 
of P 22 {t) for the cases : (i) the two atomic orbitals are ini¬ 
tially unoccupied (black long-dashed line) and (ii) both of 
the atomic orbitals are initially occupied (the green dash-dot 
line). In plot (b), the molecule starts with the one-particle 
state |1). The blue solid line with marked white circles shows 
the subsequent time evolution of pn (t) while the red solid line 
is for P 22 (t). In this figure, we use fl = SP. 


Using the above solution, we have studied the transient 
dynamics of artificial molecules triggered by different ini¬ 
tial states. We show that some initial states lead to more 
efficient occupations of the particular target states than 
the other initial states do. This also reflects the different 
pathways the molecule transverses in the presence of the 
dissipative contact reservoirs. By comparing the tran¬ 
sient probabilities of reaching one particular state from 
different initial states, one can distinguish the internal 
and the external transition pathways. 

State preparations and read outs of engineered Hamil¬ 
tonian that can mimic molecules are feasible using cold 
atoms with magneto- and optical controlsi^r.^. The 
readings of the state of the DQD molecule are often 
via separate mesoscopic electrometers 41 or alter¬ 
natively by integrating the DQD with radio-frequency 
resonant circuitsiH The quantum state properties inves¬ 
tigated with the example of artificial molecules here are 
readily to be examined in experiments. 

The general formulation presented in this article is 
applicable to time-dependent Hamiltonian. The time- 
dependence, besides arising from external fields, can also 
come from two-body interactions approximated by suit¬ 
able mean-field treatments. The general analytical so¬ 
lution to exact fermion master equations should have 
more other potential applications to various realizations 
of nanoscale quantum devices. An understanding of the 
general behaviours of the quantum states is also helpful 
to the development of quantum controls. 


In conclusion, we have obtained the general analyt¬ 
ical solution to the exact fermion master equation pre¬ 
sented in Refs. Such solution is given by Eq. da¬ 

ft connects the initial state of the principal system to 
its state at later times. It is completely in terms of 
single-particle propagating and correlating Green func¬ 
tions given by Eqs. (iTal) and dzEi, respectively. Direct 
and complete characterization of the quantum states in 
real time for a class of nanoelectronic systems and their 
cold-atom analogs is then made feasible. 
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